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Abstract 

Background: The wild barley Hordeum chilense fulfills some requirements for being a useful tool to investigate 
the endosperm yellow pigment content (YPC) in the Triticeae including its diploid constitution, the availability of 
genetic resources (addition and deletion stocks and a high density genetic map) and, especially, its high seed YPC 
not silenced in tritordeums (amphiploids derived from H. chilense and wheat). Thus, the aim of this work was to 
test the utility of the H. chilense genome for investigating the YPC in the Triticeae. 

Results: Twelve genes related to endosperm carotenoid content and/or YPC in grasses (Dxr, Hdr [synonym ispH], 
Ggppsl, Psy2, Psy3, Pds, Ids, e-Lcy, b-Lcy, Hyd3, Ccdl and Ppol) were identified, and mapped in H. chilense using rice 
genes to identify orthologs from barley, wheat, sorghum and maize. Macrocolinearity studies revealed that gene 
positions were in agreement in H. vulgare and H. chilense. Additionally, three main regions associated with YPC 
were identified in chromosomes 2H ch , 3H ch and 7H ch in H. chilense, the former being the most significant one. 

Conclusions: The results obtained are consistent with previous findings in wheat and suggest that Ggppsl, Zds and 
Hyd3 on chromosome 2H ch may be considered candidate genes in wheat for further studies in YPC improvement. 
Considering the syntenic location of carotenoid genes in H. chilense, we have concluded that the H ch genome may 
constitute a valuable tool for YPC studies in the Triticeae. 
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Background 

Hordeum chilense Roem. et Schultz. is a diploid wild 
barley included in the section Anisolepis Nevski, native 
to Chile and Argentina. This species has been used as a 
new source of cytoplasmic male sterility for hybrid 
wheat production [1-3] and it shows other potentially 
useful traits for wheat breeding, including immunity to 
Septoria tritici blotch and high seed yellow pigment con- 
tent (YPC) [4-6]. The compatibility of the H. chilense 
genome with Triticum spp. is high, and fertile amphi- 
ploids, named tritordeums (xTritordeum Ascherson 
et Graebner), have been developed from the hybrid 
between H. chilense and durum or common wheat [7,8]. 
Tritordeums are subjected to a breeding program at 
IAS-CSIC conducted by Prof. A. Martin with the goal of 
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obtaining a new crop. Additionally, they are also useful 
as a genetic bridge to transfer H. chilense genes to both 
common and durum wheat [9]. 

One of the most valuable traits of this new cereal is 
the high endosperm YPC given by H. chilense genes. 
The YPC, mainly caused by carotenoids, is also very 
important in durum wheat due to its relation to semo- 
lina quality [10-13]. In addition, the consumption of 
carotenoid-rich foods has been associated with a reduced 
risk of developing certain types of diseases [14,15]. 
Besides, carotenoids with unsubstituted (3-ionone rings, 
such as |3-carotene, have provitamin A activity [16] while 
other carotenoids such as lutein may prevent age-related 
macular degradation [17,18]. Accordingly, seed caroten- 
oid content is receiving more attention with the goal of 
developing new functional foods. 

The carotenoid biosynthesis pathway is well known in 
plants [16,19-22]. Candidate gene studies in maize have 
shown the role of several carotenoid-related genes in 
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endosperm carotenoid content using quantitative trait 
locus (QTL) and association studies. Indeed, several 
genes including Psyl [23-25]; Zds [25]; e-Lcy [26]; Ziso 
[27] and Hyd3 [28] have been related to endosperm 
carotenoid content. Genes from both upstream and 
downstream pathways have also been investigated at the 
transcriptomic level in maize [29]. 

Despite the importance of YPC in durum wheat and 
tritordeums and the positive correlation between YPC 
and carotenoid content, candidate gene studies for YPC 
in both species are mostly based on Phytoene synthase 1 
gene (Psyl). PSY1 catalyzes the first step of the caroten- 
oid biosynthetic pathway and it is considered a limiting 
factor for carotenoid production [18,30]. Since the first 
reports about the location of Psyl to chromosomes 7 A 
and 7B in wheat [31,32], a considerable amount of work 
has demonstrated the role of this gene in durum wheat 
YPC [33-39]. In addition, the lycopene epsilon cyclase 
gene (e-Lcy) has been co-localized with a QTL for YPC 
in wheat [35]. Nevertheless the role of other candidate 
genes for YPC in durum wheat remains largely unex- 
plored. Given the high level of macrocolinearity among 
grasses, the knowledge gained in any species may be 
useful in others. This is more evident within the Triti- 
ceae due to the close relationship among these species. 
Indeed, Triticeae species are represented as a single 
genome in comparative studies [40,41]. 

In this context H. chilense may be a useful species to 
investigate the role of carotenoid-related genes in the 
YPC in the Triticeae. First, it exhibits a high YPC which 
is not silenced in wheat background, as evidenced in tri- 
tordeums [10,42,43]. Second, the diploid constitution of 
this species and the existence of H. chilense-wheat 
chromosome addition lines [44] and other deletion 
stocks [45,46] makes the identification of candidate 
genes more accessible. Third, a high density genetic map 
has been recently developed using H. chilense-derived 
DArT markers [47] which allows efficient candidate gene 
mapping and comparative studies. 

The objective of this work is to determine whether 
H. chilense could be useful for investigating the endo- 
sperm YPC in the Triticeae. For this purpose we aimed 
to identify and locate a set of carotenoid-related genes 
previously studied in other grasses and to determine the 
macrocolinearity between H. chilense and barley for 
these genes. In the second place, we aimed to identify 
the main regions for YPC in H. chilense and to discuss 
the results obtained in relation to wheat. 

Methods 

Plant material 

Hordeum chilense accession lines HI and H7 from 
the collection of the Institute de Agricultura Sosten- 
ible (IAS-CSIC, Cordoba) were used for polymorphism 



detection in the selected genes. DNA of 92 F 7 -Recom- 
binant Inbred Lines (RILs) derived from the cross HI x 
H7 was used for genetic mapping of the gene-based 
markers designed. These lines were sown in small pods 
with five seeds each. After initial growing, RILs were 
transplanted to field conditions using a randomized plot 
design with two replications. 

For physical mapping, the following common wheat 
cv. Chinese Spring (CS)-//. chilense addition lines were 
used: with complete chromosomes lH ch , 4H ch , 5H ch , 
6H ch or 7H ch (named CS MA lH ch -l H ch S, CS DA4H ch , 
CS DA5H ch , CS DA6H ch and CS DA7H ch , respectively, 
where MA refers to monosomic addition and DA to 
disomic addition), and (CS)-//. chilense ditelosomic 
addition lines CS DAlH ch S, CS DA2H ch S, CS DA5H ch L, 
CS DA6H ch S, CS DA7H ch a, CS DA7H ch [3. 

Leaf tissue was harvested at tillering stage, frozen in 
liquid nitrogen, and stored at -80°C. Genomic DNA was 
extracted using the CTAB method according to [48]. 

Primer design 

The twelve genes were selected from the main steps of 
the carotenoid biosynthetic pathway but also from up- 
stream and downstream pathways, considering previous 
results obtained in other grasses, mainly maize. Rice 
genes were used as a query to identify orthologous genes 
in maize, sorghum, wheat and/or barley (see Additional 
file 1). Nucleotide sequences were aligned to design 
primers in the conserved regions of each gene. Psyl 
sequence of H. chilense was available from previous 
works [31,49]. Sequences were aligned using Edialign 
program (http://emboss.sourceforge.net/index.html) and 
edited using GeneDoc software (http://www.psc.edu/ 
biomed/genedoc). Sequence identity searches were per- 
formed at the NCBI (http://www.ncbi.nlm.nih.gov) using 
BLAST. Primer pairs were designed using Primer3plus 
software [50] on exonic regions flanking at least one in- 
tron when possible. 

PCR amplification and polymorphism detection 

PCR amplifications were performed in 25 ul reactions 
consisting of 0.625 units of DNA polymerase (Biotools 
B&M Labs, Madrid, Spain), lx PCR buffer, 1.6 mM 
MgCl 2 , 320 mM dNTPs (Promega, Madison, WI, USA), 
0.6 mM of each primer and 50 ng of genomic DNA. 
PCR were carried out as follows: 5 min at 94°C, 35 cycles 
of 30 s at 94°C, 30 s at 58°C (54°C for Zds) and 1 min 
at 72°C followed by 7 min at 72°C. For allele-specific 
amplifications with tetra-primer PCR, the number of 
cycles was increased to 40 and the annealing tem- 
perature was set at 60°C. Amplification and digestion of 
Psyl gene were carried out as described by [31]. PCR 
products were resolved on 2% agarose gels, stained with 
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ethidium bromide and visualized under UV light. All 
primers used in this work are summarized in Table 1. 

Nucleotide sequences were obtained by direct sequen- 
cing or by cloning if more than one product was 
observed in the amplifications. For cloning, pGEMT- 
Easy vector (Promega, Madison, WI) was used to trans- 
form competent Escherichia coli (DH5a) cells. Plasmids 
were isolated and purified using Illustra plasmid Prep 
Mini Spin Kit (GE Healthcare, UK Ltd, UK). 

Digestion reactions with restriction enzymes Hhal, 
Hinfl, MboII, Nhel and TasI (Takara Bio Inc., Japan) 
were carried out following manufacturer's instructions 
using 10 ul of PCR in a final volume of 20 ul. 

Identification of genomic regions for endosperm yellow 
pigment content (YPC) 

YPC was determined according to the AACC method 
14-50 [51] in 71 F 7 -RILs, which yielded a sufficient 
amount of seed. Briefly, approximately 150 mg of 
flour was used to evaluate YPC. Pigment extraction 
was performed using water-saturated n-butanol in a 
1:10 (w/v) ratio. Samples were centrifuged at 14,000 g 
for 10 min and the pigment was measured by absorb- 
ance reading at 435.8 nm. All determinations were per- 
formed in duplicate. 

A DArT-based H. chilense map [47] constituted the 
basis for mapping. A framework map composed by mar- 
kers spaced 3-5 cM was considered to determine trait- 
marker associations and QTL mapping since the use of 
very proximal markers results in highly similar LOD 
scores while it signifies a huge increase in the number of 
calculations. Gene-based markers were integrated in the 
H. chilense map using JoinMap 4.0 mapping software 



[52]. After linkage analysis, the position of each candi- 
date gene within a chromosome was determined using 
the fixed order option, the Monte Carlo Maximum Like- 
lihood Mapping algorithm (ML-mapping) module and 
the Kosambi mapping function. 

Associations between markers and YPC were deter- 
mined using the Kruskal-Wallis (KW) test and MapQTL 
4.0 [53]. Regions showing at least five markers associated 
with YPC at p<0.01 were considered as YPC regions but 
further analyses were performed to identify the most 
significant ones. Marker cofactors were automatically 
selected using the automatic cofactor selection (ACS) 
tool at p<0.05. For both Interval Mapping (IM) and Mul- 
tiple QTL Mapping (MQM), the genome-wide signifi- 
cance threshold was estimated using the permutation 
test with the following parameters: the walk speed was 
set to 3 cM; the number of permutations was set to 
500 and the significance level to 0.05. The QTL position 
was estimated at the peak LOD with a 2-LOD support 
interval [54]. The map figure included in the manuscript 
was produced using MapChart software [55]. 

Results 

Gene-based marker design, amplification and mapping 

A total of twelve genes not previously mapped in 
H. chilense were selected for this work: 1-deoxy-D- 
xylulose 5-phosphate reductoisomerase (Dxr) and 4- 
hydroxy-3-methylbut-2-enyl diphosphate reductase (Hdr 
[synonym ispH]) from the MEP pathway; geranyl geranyl 
pyrophosphate synthase (Ggppsl) for geranylgeranyl 
diphosphate synthesis; phytoene synthase 2 (Psy2), phy- 
toene synthase 3 {PsyS), phytoene dehydrogenase (Pds), 
zete-carotene desaturase {Zds), lycopene epsilon cyclase 



Table 1 Primers designed for amplification and mapping of the candidate genes 



Primer 


Sequence 5'-3' 


Primer 


Sequence 5'-3' 


GGPPS1-F1 


GCCAGCGTCGACTCCTAC 


PSY3-F 


TCGACGAGCTGTACGTCTACTG 


GGPPS1-R1 


GGAACAGCAAGCCAATTGAT 


PSY3-R 


CTCGATCTCGTCCAGGATCT 


GGPPS1-spF 


CGTATGGCTTFCTAAGAAGTG 


HYD3-F3 


CCTCCGTGTACTAGCGCTTG 


GGPPS1-spR 


GGAGATACCTATGCAAATCAT 


HYD3-R2 


CGAAGTFGTGCATGTGGTGT 


DXR-F1 


GAGCATGGGGAAGAAAATCA 


ZDS-F2 


TGTCCCAGGGATGAAAAGAC 


DXR-R1 


GGTGACCTCGGAGCAATAGA 


ZDS-R2 


AGCTTGGATCAGGGAACCTT 


HDR-F2 


GGCATTGCAAATCAAACAAC 


LBG-F2 


CAAGGTGAAGTCCACCATGA 


HDR-R2 


TTTCCTGGTGCAATGCTTTG 


LBC-R2 


TGTCTAGAAACCGCACGATG 


LEC-F2 


GTGGACAATATTTGCCTGGAA 


LBC-spF 


GCGACTCCCACCTCCCT 


LEC-R2 


GGAGGTGTGTGAGGAGGTFC 


LBC-spR 


GATGGCGGACCGGTCG 


PDS-F1 


TAGAGGTCGTGATTGCTGGT 


CCD-F1 


TTGATCCTACAAAGAAAGGTGGT 


PDS-R1 


GGGAAATCAAAACGGCTGTA 


CCD-R1 


CAGCTCATTCGCGAAGTTCT 


PSY2-F 


ATTGCTCCGGACTCAAAGG 


CCD-spR 


GAAGTTCTCGAGCTTGTCGG 


PSY2-R 


TTGTAGTCGTTGGCGTGGAT 


PP01-F1 


AGCTTGGAGGAGCAGTGG 


PSY2-spR 


GAGTGAACAATGCTTGAATGA 


PP01-R1 


GTGGTGCGCGAAGAAGAT 
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{e-Lcy), lycopene beta cyclase (b-Lcy) and beta-carotene 
hydroxylase 3 (HydS) from the carotenoid biosynthetic 
pathway; the carotenoid cleavage dioxygenase 1 gene 
(Ccdl) mediating the degradation of carotenoids to apoc- 
arotenoids and polyphenol oxidase 1 gene {Ppol) impli- 
cated in plant tissue enzymatic browning. 

Rice genes (www.gramene.org) were used as a query to 
identify orthologs from barley, wheat, sorghum and 
maize. Exon-intron boundaries and highly conserved 
regions were recognized by the alignment of each 
sequence set. Primers were designed on conserved exo- 
nic regions to amplify at least one intron and maximize 
polymorphism detection. Amplifications with the pri- 
mers designed (Table 1) were carried out using the par- 
ental lines HI and H7 as template. PCR products were 
sequenced to confirm their identity and examined for 
polymorphisms by comparing both parental lines. 

All the candidate genes were successfully amplified 
and sequenced in H. chileme [GenBank: JQ922078- 
JQ922105]. Differences between HI and H7 sequences 
were found for all the candidate genes, these being in 
most cases Single Nucleotide Polymorphisms (SNP) 
(Table 2). These SNPs were used to develop either 
Cleaved Amplified Polymorphism (CAP) markers for 
Zds, HydS, Pds, Hdr, Psy3 and e-Lcy or to design allele- 
specific primers. Specific amplifications were obtained 
by using tetra-primer PCR (to amplify both HI and H7 
alleles) for Ggppsl (GGPPS1-F1/ GGPPS1-R1/ GGPPS1- 
spF/ GGPPSl-spR) and b-Lcy (LBC-F2/ LBC-R2/ LBC- 
spF/ LBC-spR). Alternatively, specific reverse primers 
CCD-spR and PSY2-spR were designed to use in com- 
bination with CCD-F1 and PSY2-F, respectively, to 
obtain PCR amplification exclusively in one of the 



parental alleles. Finally, length polymorphisms were 
detected for Dxr (14 bp insertion in H7 sequence com- 
pared to HI) and Ppol. Regarding Ppol, two sequences 
were amplified in HI with primer pair PPOl-Fl/PPOl- 
Rl differing in two SNPs (JQ922103 and JQ922104), 
both of them having in common a deletion of 23 bp 
compared to H7 (JQ922105). This was the length poly- 
morphism scored in the mapping population. 

All gene-based markers were located in the H. chilense 
map (Figure 1) in the following chromosomes: 2H ch 
{Ggppsl, Zds, Hyd3 and Ppol); 3H ch {Dxr and e-Lcy); 
4H ch {Hdr and Pds); 5H ch {Ccdl, Psy2 and Psy3); 6H ch 
{b-Lcy). The position (cM) of each candidate gene in the 
map is shown in Table 3. Both Ggppsl and Zds were 
tightly linked and mapped near the centromere which 
was estimated in a 11 cM window [47] while Hyd3 
and Ppol were mapped to the long arm of chromo- 
some 2H ch at 142.38 cM and 161.56 cM respectively. 
The e-Lcy gene was located in the vicinity of the 
chromosome 3H ch centromere while Dxr was mapped to 
the distal part of the 3H ch S arm (Figure 1). Hdr and Pds 
genes were mapped to 4H ch S and 4H ch L arms respect- 
ively. Ccdl (distal) and Psy2 (39.06 cM) were mapped to 
chromosome 5H ch S while Psy3 was mapped to the long 
arm of this chromosome (133.75 cM). Finally, b-Lcy was 
mapped to the long arm of chromosome 6H ch (111.2 cM) 
since the centromere is estimated at 84.0 cM [47]. Also, 
Psyl is located on chromosome arm 7H ch a as shown in 
previous studies using physical mapping [31,47]. 

Physical mapping by using the series of addition lines 
of H. chilense in common wheat was also possible 
for genes Hdr, Ggppsl, Psy2, Zds and Ccdl. Physical 
location could not be confirmed for genes located on 



Table 2 Polymorphism between HI and H7 lines used for mapping 



Gene 


Accession numbers 




Polymorphism detection 




HI 


H7 


Genetic mapping 


Physical mapping 


Dxr 


JQ922078 


JQ922079 


Allele size 




Hdr 


JQ922096 


JQ922097 


Mboll digestion 


Mboll digestion 


Ggpps 1 


JQ922098 


JQ922099 


Allele specific primers 


Allele specific primers 


Psy2 


JQ922084 


JQ922085 


Allele specific primers 


Allele specific primers 


Psy3 


JQ922082 


JQ922083 


Hhal digestion 




Pds 


JQ922092 


JQ922093 


Hhal digestion 




Zds 


JQ922090 


JQ922091 


Hinfl digestion 


Tasl digestion 


e-Lcy 


JQ922094 


JQ922095 


Hhal digestion 




b-Lcy 


JQ922088 


JQ922089 


Allele specific primers 




Hyd3 


JQ922080 


JQ922081 


Nhei digestion 




Ccdl 


JQ922086 


JQ922087 


Allele specific primers 


Allele specific primers 


Ppol 


JQ922103 


JQ922105 


Allele size 






JQ922104 









Polymorphism between H1 and H7 lines of H. chilense used 
number for each sequence is given. Psyl is not included as 



for genetic and physical mapping of the candidate genes selected in this work. GeneBank accession 
sequencing and mapping have been reported elsewhere. 
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1H 



ch 



2H 



ch 



3H 



ch 



4H 



ch 



5H 



ch 



6H 



ch 



7H 



ch 



- bPb-791825 
-bPt-790798/1Hch 
-bPt-34026/1Hch 
-bPt-801126/2-3-1Hch 

- bPt-790498 

- bPb-20089 
J<N bPb-803863 

^ bPt-788790/H7 
- ^bPt-790562/1Hch 



■bPb-801070 

- bPt-787804/H7 

- Cos250 
-bPt-801106n"21 
-bPt-39113/1Hch 
-bPb-787817 

- Cos252 

- bPt-789995 
-S208.1 

- bPb-43060 
-bPt-791 247/1 Hch 
-bPb-801860 

- bPb-23208 

- bPb-803622 
-bPb-801220 
-bPb-802160 

- bPt-789909/H7 
-bPt-801 928/1 Hch 
-bPb-802819 



-bPt-792286/1Hch 
-bPt-804165 
-bPt-790546/1Hch 
-bPb-801304 
- bPt-803708 



^bPt-801 084/1 Hch 



r bPb-788872 
^bPt-791659/2Hch 
^bPt-791350/2Hch 
^ bPb-790087 
^ bPb-803693 
^ bPt-792059 
-bPt-801503/2Hch 

- bPt-803786 

- bPt-804635 

- bPb-802898 
^167F-IAS-pHc2-1 
^ bPb-803080 



^ bPt-803049 
^bPb-801877 
^ bPt-801598 



bPI-787052 

- S208.4 
•-bPt-789148/H7 
-bPt-801612 

- Hyd3 
-bPb-801888 
^ bPt-787983 



- bPt-27784/2-3Hcb 

- bPt-800966 

- bPb-804306 

- bPt-790051 

- bPt-804393 

- bPt-786845/H7 
--bPb-801217 

- bPt-791531 

- bPt-788067 

- bPt-788984/H7 

- bPt-789975/H7 
-bPb-801175 

- bPt-790317 
bPb-804457 



i 



-bPb- 14883 

- bPI-787861 
-bPb-804152 
•-bPt-801920/2-3Hc 
^bPt-14662 
-bPb-791376 
-bPb-787128 
-bPt-801024/2-3Hc 
-bPt-787128 
-bPb-801449 

- bPb-800832 
r bPb-802358 

bPt-790363 
bPt-793041/2-3Hc 
bPt-787249 
bPt-792785/2-3Hch 
bPt-792 3477T21 

_NbPt-792931/2-3Hc 
^ bPb-803446 

_-bPb-23391 

_^bPt-792986/2-3Hch 

_^bPb-802564 

_^bPt-0959 

--bPt-22320 

_^S197 

bPb-804555 

bPt-790460/2-3Hc 

z-bP1-801075/2-3Hc 
_^bPb-802382 
_ S bPt-788034 
-^bPb-803232 

-S bPt-791546 

_^bPb-791293 

^bPb-801924 
---bPb-786929 

^bPb-802154 
-^bPb-5245/4H-xx 
~^bP1-801466/2-3Hc 
_^bPb-801864 
= — bPt-48561 
" bPb-787055 



- bPb-804484 

- bPt-788698 

- bPt-787644 
-bPt-800818/4Hch 

- bPt-22046 
-S193 

- bPt-790155 
bPb-803647 

- bPt-789229 

■H 

- bPI-788912 
-bPb-801750 
•-bPt-787419/4Hch 
^ bPt-789551/H7 

bPt-789475/H7 
Cos286 

bPt-792793/4Hch 



SB 

\V bPt-7 



bPt-792250/4Hch 
bPb-790810 
bPt-803175 
bPb-802140 
bPt-789418/H7 
bPb-804636 
bPt-49174 
bPt-790739 
bPb-804518 
bPt-38224 
bPt-35303/H7 



^bPb-801384 
f bPb-802097 
_<CbPb-49313 
->bPM814 
J,Psy2 

J, bPb-803321 

- / bPb-788862 
= i-bPt-790333 

■ , bPt-37496/5Hch 
■7bPt-786619 
■/-bPt-801758 
-ybPt-789395 
bPt-790102 

- , hPI-7Q144D 



-bPt-787198 

- bPt-787587 

- Cos293 



bPt-789102 
bPt-786654 
V-bPb-788614 
\ v bPb-802202 
bPt-789006 
\ v bPt-45020 
\ v bPt-803799 
A bPt-790576/5Hch 
'O-bPt-12411 
'O-bPb-791887 
YbPt-790677 
^bPt-46182 
bPt-791836 
bPt-787206/5Hch 
bPt-788057 



bPt-789851/H7 
- Bmac31 6 
bPf-802289/6Hch 
bPf-787342 
bPt-47201 
bPb-B03627 
bPf-802544/6Hch 
bPb-788586 
bPf-790540/6Hch 

■ bPb^5053 

■ bPt-791425 

■ bPt-791407 

■ bPt-791866 

■ TAFAD 

bPt-47517/6Hch 
bPf-792043/6Hch 
-bPt-33 143/6Hch 

I, S183.1 
l=V bPt-788758 
/■ bPb-803976 
, bPt-45095 
f bPb-802956 
^-bPt-788106/H7 
^bPb-791305 
r bPb-8Q0982 

- bPt-787666 

- bPb-789364 
-bPt-38718 
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Figure 1 Marker-YPC association in the Kruskal-Wallis (KW) analysis. To determine marker-YPC associations, a framework map composed by 
markers spaced around 3-5 cM was considered from the map developed by [47]. The estimated centromere position is shown as a coloured 
segment within each chromosome. In chromosomes 6H ch and 7H ch centromeres were located in a single position since markers belonging to 
different chromosome arms mapped at the same position. For chromosomes 6H ch and 7H ch , the estimated centromere locations are shown as 
black triangles. The position of the candidate genes mapped in this work is highlighted using a coloured box. The significance of the KW analysis 
is shown: - Not significant; *0.1; **0.05; ***0.01; ****0.005; *****0.001; ******0.0005; *******0.0001 . Three locations with at least five markers 
associated with YPC were considered as YPC regions (Regions 1 to 3). The most significant zone within these regions was determined using 
MQM analysis and permutation test (identified as YPC at 2H ch L) and its position estimated at the peak LOD with a 2-LOD support interval. 



chromosomes 2H C L and 3H C (since there are no 
addition lines for these chromosomes) or when poly- 
morphism with wheat was not detected with the 
designed markers (b-Lcy y Psy3, Pds and Ppol). 

To study the macrocolinearity between H. chilense 
and barley for these candidate genes, the Poaceae ortho- 
logous group corresponding to each candidate gene 
was identified from http://rice.plantbiology.msu.edu/ 
annotation_pseudo_pog.shtml using the rice gene as 
query (Table 3). Orthologous genes in Brachypodium 
and Sorghum were identified by this means. The RAP 
locus identifier was retrieved using the ID Converter 
tool (http://rapdb.dna.affrc.go.jp/tools/converter). Rice, 
barley and sorghum orthologs were located in the Barley 
Genome Zipper [56], which allowed the determination 
of their relative positions in the barley map [57] in all 
cases except for Zds which was not found in the Barley 
Genome Zipper (Table 3). The Psyl, previously mapped 
in H. chilense, was also located in barley. 

Determination of YPC regions in H. chilense 

The YPC content was determined in the HlxH7 RIL 
mapping population. Kruskall-Wallis (KW) association 
test revealed three main chromosome regions related to 
YPC (Figure 1). The largest is located on chromosome 



2H C where the majority of markers showed an associ- 
ation with YPC (Figure 1). Smaller regions were detected 
on chromosomes 3H ch and 7H ch . Candidate genes 
Ggppsl, Zds, Hyd3, Ppol, e-Lcy and Psyl showed a 
significant association with YPC (Figure 1). The MQM 
analysis resulted in a significant QTL at the distal part 
of chromosome 2H ch between 207 and 229 cM O?<0.05) 
where no carotenoid-related genes were mapped in 
this work. 

Discussion 

Macrocolinearity between H. chilense and H. vulgare 

Comparative analysis of the twelve candidate genes 
positions in H. chilense and barley using the Barley 
Genome Zipper [56] revealed a good agreement between 
both species for the candidate genes selected (Table 3). 
Indeed, chromosome location and relative positions of 
genes within chromosomes were conserved in all cases. 
Regarding Zds, its position could not be established in 
barley since the orthologous genes from rice, sorghum 
or Brachypodium were not found in the Barley Genome 
Zipper [56]. However, the close linkage between Zds 
and Ggppsl in H. chilense suggests that the two genes 
are located in equivalent positions in barley. This is 
supported by the high-resolution comparative analysis 
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Table 3 Macrocolinearity between H. chilense and H. vulgare for the candidate genes selected 
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Identifier 1 


Rice LOC 2 


Rice RAP 3 
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2 
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e-Lcy 


POG. 


ORTHOMCL17555 
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Sb03g026020 


Bradi2g41890 


3 


55.57 


3H ch C-L 


1 00.47 


Dxr 


POG_ 


ORTHOMCL8746 


LOC_Os01y01710 


Os01g01 06900 


Sb03g008650 


Bradi2g00650 


3 


0 


3H ch S 


7.45 


Pds 


POG. 


ORTHOMCL8451 


LOC_Os03g08570 


Os03g01 84000 


Sb06g030030 


Bradi1g72400 


A 


89.39 


4H ch L 


83.39 


Hdr 


POG. 


ORTHOMCL5875 
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Sb01g009140 
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4 
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4H ch S 


47.49 
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LOC_Os09g38320 
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5 
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6 
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6H ch L 
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Psy1 
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C) Poaceae orthologous group corresponding to each candidate gene; ( 2 ) MSU rice database identifiers of rice genes used as query to identify the Poaceae 
orthologous group; ( 3 ) Rice, Sorghum and Brachypodium genes included in each Poaceae orthologous group. For nomenclature of rice genes (LOC vs RAP) see 
(http://rapdblegacy.dna.affrcgo.jp/gene/nomenclature.html); ( 4 ) Position of the candidate genes in barley according to the Barley Genome Zipper. ( 5 ) Position 
of the candidate genes in hi. chilense map. Genes in bold are shown in the Barley Genome Zipper [56]; nd: not determined. 



between barley and rice performed by [56]. Briefly, 
the Barley Genome Zipper establishes the putative lin- 
ear order of 21,766 barley genes by integrating gene 
indices of rice, sorghum and Brachypodium [56] and 
allows the establishment of the synteny relations 
between barley and rice chromosomes based on the 
position of the putative orthologs in both species. 
Chromosome 2H in barley corresponds to chromo- 
some 4R in rice, except for a region in 2HS arm which 
is equivalent to 7R. In barley the region near the 2H 
centromere shows significant matches to both 7RS 
(where Zds is located in rice) and 7RL (where Ggppsl 
is located in rice). This is in agreement with H. chi- 
lense results and suggests that Zds would be located in 
an equivalent position in barley. 

An apparent break in the synteny between H. chilense 
and H. vulgare happens for Psyl which is located in 7HL 
in barley and in 7H ch oc in H. chilense [31]. However, 
genetic mapping could not establish the relative position 
of Psyl in 7H ch chromosome since it remained in a short 
linkage group not linked to 7H ch chromosome. Thus, its 
putative location in 7H ch a is based upon physical map- 
ping using a telosomic addition line. We have recently 
observed a reorganization of the H. chilense chromatin 
in wheat background between different chromosomes 
(A. Martin, personal communication). Thus, a transloca- 
tion involving Psyl and 7H ch S during the development 
of the addition line cannot be totally excluded and fur- 
ther work is required to confirm the position of Psyl in 
the chromosome 7H ch . 



YPC regions: comparative analysis between H. chilense 
and wheat 

The association study for YPC allowed the identification 
of three main regions in chromosomes 2H ch , 3H ch and 
7H ch for YPC variation in H. chilense as revealed by a 
KW association test. Despite a single QTL at the distal 
part of chromosome 2H ch L being above the significance 
threshold established, all three regions were considered 
significant since the population size critically affects the 
power of QTL detection and the precision of QTL esti- 
mates [58]. For instance, QTLs with small effects for 
stripe rust resistance in barley were detected only by 
increasing the population size [59]. The significance 
threshold must thus be chosen according to the popula- 
tion size and the trait of interest [58]. Therefore, we 
relied on KW-defined regions and considered the QTL 
analyses to conclude that the distal part of chromosome 
2H ch was the most significant one. Seven of the candi- 
date genes were not associated with YPC in this work. 
This may indicate that they do not contribute to YPC 
variation in this population. Alternatively, these genes 
might have a minor effect which cannot be detected with 
this population size. 

The region located in chromosome 3H ch includes the 
e-Lcy gene, whose role in YPC has been demonstrated in 
hexaploid wheat and maize. In wheat, e-Lcy co-localized 
with a QTL for lutein accumulation and a specific muta- 
tion causing lutein content variations was identified [35]. 
This QTL was located in the proximal 3BS region this 
being the position of this gene in hexaploid wheat in 
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agreement with the results obtained in H. chilense. Simi- 
larly, it has been shown that a QTL for carotenoid accu- 
mulation in maize chromosome 8 co-locates with e-Lcy 
gene [25,26]. 

The distal part of chromosome 7H ch L is also related 
to YPC in H. chilense. The importance of this region 
has been clearly demonstrated in wheat, as a result of 
variations in Psyl genes [32,33,38]. This has promoted 
extensive work towards the identification of Psyl-A and 
Psyl-B alleles, the determination of their role in YPC 
and the development of functional markers for selection 
[34,60,61]. Furthermore, there is evidence of an add- 
itional locus for YPC in the distal part of chromosome 
7AL in durum wheat [38]. Given that Psyl is apparently 
located in the 7H ch a arm, this might indicate that an 
additional locus is found at 7H ch L, in agreement with 
durum wheat results. Nevertheless, since genetic map- 
ping did not resolve Psyl position within chromosome 
7H ch , we cannot totally discard the occurrence of a 
translocation 7H ch L/7H ch S during the development of 
the wheat///, chilense chromosome addition line. Thus, 
macrocolinearity cannot be established for this gene at 
present and the YPC region on chromosome 7H ch L 
might possibly be related to Psyl variation. Indeed, Psyl 
showed a low association with YPC (p<0.l). 

Finally, there is a region on chromosome 2H ch related 
to YPC which extends over the 80% of the chromosomal 
length. Four candidate genes, Ggppsl, Zds, Hyd3 and 
Ppol, are mapped within this interval. The most signifi- 
cant zone within this region is located in the distal 
part of chromosome 2H ch L (Figure 1). However, the 
complete KW-defined region must be considered as a 
target since the population size is probably hampering 
the determination of other QTLs in this chromosome. 
Indeed, in a preliminary map constructed using a 
HlxH7 F 2 population, which did not cover the distal part 
of chromosome 2H ch , a QTL for YPC was detected in the 
chromosome 2H ch in the vicinity of the marker IAS- 
pHc2-l [62] which is now mapped at 48 cM (Figure 1). 
Additionally, the location of four carotenoid-related genes 
in this area further supports the importance of this region 
for YPC. 

The relevance of chromosome 2A for YPC and/or 
carotenoid content has been previously reported in 
wheat. Most works have focused on the role of polyphe- 
nol oxidases since these enzymes cause the darkening of 
flour, pasta and noodle products [63]. Duplicated poly- 
phenol oxidase genes {Ppol and Ppo2) have been identi- 
fied and mapped clustered in barley 2H chromosome 
[64]. Although only Ppol was mapped in H. chilense, a 
tightly linked Ppo2 may be expected. In any case, the 
association of Ppol and YPC in H. chilense seems to 
indicate the occurrence of darkening processes in this 
species. Similarly, the role of Ppo-Al and Ppo-Dl has 



been extensively studied in wheat including mapping 
studies [65,66], the development of functional markers 
[67-69] and the cloning and analysis of these genes 
[70,71]. Recently, duplicated Ppo genes (Ppo-A2, Ppo-B2 
and Ppo-D2) have been identified [72] and mapped to 
chromosomes 2A, 2B and 2D [73]. Thus the results 
obtained in H. chilense are in agreement with previous 
findings in wheat. 

Our results also suggest the potential importance of 
the centromeric region on chromosome 2H ch for YPC. 
QTL and association mapping studies have identified 
regions for YPC or carotenoid content in chromosome 
2A [32,33,74]. Recently, a QTL on the centromeric 
region of chromosome 2 A has been reported [33]. How- 
ever, it could not be associated with YPC per se since the 
markers within the confidence interval were related to 
kernel weight variation. Thus, it was proposed that this 
QTL may be associated with yellow pigment concentra- 
tion due to a lower carbohydrate content in smaller 
grains or it could be linked to a QTL for kernel weight 
[33]. Since there is no significant association between 
YPC and kernel weight in H. chilense [42], this region 
does not seem to be related to carotenoid concentration. 
Moreover, previous findings in durum wheat also sup- 
port the existence of a QTL for YPC in this region 
[32,74]. Both studies showed the association between 
Xgwm49S and YPC using a mapping population derived 
from the cross W9262-260D3 x Kofa [32] and a col- 
lection of 93 diverse accessions of durum wheat [74]. 
This marker is located in the arm 2AS (C-2AS5-078 bin) 
at 108.5 cM [75] (see Group 2 data file at http://wheat. 
pw.usda.gov/ ggpages/SSRclub/GeneticPhysical/) between 
markers Xgwm372 (C-2AS5-078 bin, at 90.8 cM) and 
Xgwm249 (C-2AL-1-0.85 bin, at 111.9 cM). Both 
Xgwm372 and Xgwm249 are within the confidence inter- 
val for the QTL on chromosome 2A reported by [33]. 
Hence, all three works identify the same region around 
the centromere of chromosome 2A, and thus, this region 
seems to be effectively related to YPC and/or carotenoid 
content in durum wheat. The results presented in this 
work also support the importance of the centromeric 
region of the chromosome 2H ch in YPC. Furthermore, 
the Zds and Ggppsl genes were mapped within this 
region, and they should therefore be considered in 
durum wheat. The potential association of Zds with YPC 
in durum wheat has already been hypothesized after 
physical mapping on chromosomes 2A and 2B [76]. 
Later studies cloned Zds gene in common wheat [77] 
and identified a QTL for YPC related to Zds-Dl vari- 
ation in chromosome 2D in common wheat [78]. Simi- 
larly, the co-localization of Zds with a QTL for 
carotenoid content in chromosome 7 of maize [25] also 
supports the potential role of this gene in grasses. Simi- 
larly, the association between Ggppsl expression and the 
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carotenoid content in maize [29] and the close linkage 
between Zds and Ggppsl in H. chilense make the latter 
an interesting candidate gene for durum wheat studies. 
Likewise, variations at Hyd3 locus are important for 
maize endosperm carotenoid content [28,79]. Given the 
association of this gene with YPC and its position in 
H. chilense map, Hyd3 should be also considered in fur- 
ther studies in durum wheat. 

The most significant region of 2H ch , as identified by 
MQM analysis, is located in the distal part of the 
chromosome 2H ch L where no candidate genes were 
located in this work. The occurrence of this association 
can be explained by several hypotheses, although they 
are all very speculative. The simplest one is that other 
carotenoid genes not considered in this study are located 
in the distal 2H ch L. A second possibility would be that 
genes for carotenoid esterification are located in this 
region. A QTL for lutein esters was identified in 
chromosome 2B [35]. It is flanked by the marker 
KsuD22 which has been located in the chromosome 2BL 
by deletion mapping [75]. Esterification is a common 
mechanism for sequestering carotenoids in plants [80]. 
Indeed, carotenoids are accumulated in specialized 
lipoprotein-sequestering structures in the chromoplasts 
[81] and the apolar compounds of some chromoplasts, 
such as fibrils, are mainly composed of esterified xantho- 
phylls [81]. H. chilense genome addition in tritordeums 
produces both high carotenoid content and an increase 
in lutein esterification [10], and, thus, the presence of 
genes for differential esterification ability in the 2H ch L 
chromosome might be considered in the future. 

The lack of association between the rest of the candi- 
date genes (Dxr, Hdr, Ccdl, b-Lcy, Pds, Psy2 and Psy3) 
and YPC in H. chilense suggests at least two possibilities: 
(1) these genes are not related with YPC variation in the 
population; (2) they have minor effects which cannot be 
detected with the current population size. With the 
results presented in this work it is not possible to dis- 
cern between these possibilities. While Psy2 and Psy3 
are not associated with seed carotenoid content in other 
grasses (reviewed by [82]) and thus, they may be dis- 
carded for further analyses, the rest of the genes might 
be still useful for durum wheat studies. For instance, 
Ccdl copy number is negatively correlated with caroten- 
oid content in maize [83]. It is mapped in the distal part 
of 5H ch S in concordance with QTLs for carotenoid con- 
tent and YPC in durum wheat [33,84]. Thus, it may be 
an interesting candidate gene for durum wheat despite 
the lack of association with YPC in this work. 

Conclusions 

The equivalent location of the carotenoid-related genes 
in H. chilense and barley shows a high level of collinear- 
ity between both species, at least for the genes studied in 



this work. Considering the high synteny among Triticeae 
species along with the high YPC content of H. chilense 
and its derived amphiploids (tritordeums), the H ch gen- 
ome constitutes a valuable tool for YPC and carotenoid 
content studies in the Triticeae. Indeed, the YPC regions 
detected in H. chilense showed good agreement with 
previous findings in wheat. Accordingly, the location of 
Zds and Ggppsl within one of these regions around the 
2H ch centromere is in concordance with QTL for YPC 
and carotenoid content in wheat, and suggests that both 
genes should be considered for durum wheat improve- 
ment. In addition, the identification and mapping of 
the candidate genes in H. chilense will enable further 
studies to investigate the genetic basis of the high YPC/ 
carotenoid content in tritordeum which could also be 
extrapolated to durum wheat. 

Additional file 

/ 

Additional file 1: Table SI. Accession numbers of the sequences used 
for alignments in the design of primers for amplification of the orthologous 
genes in H. chilense. 
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